Journal  of  Coastal  Research 


21 

3 

522-534 

West  Palm  Beach,  Florida 

May  2005 

Beach  Profile  Equilibrium  and  Patterns  of  Wave  Decay  and 
Energy  Dissipation  across  the  Surf  Zone  Elucidated  in  a 
Large-Scale  Laboratory  Experiment 

Ping  Wang+  and  Nicholas  C.  Kraus* 


Tlepartment  of  Geology 
University  of  South  Florida 
4202  E.  Fowler  Ave. 

Tampa,  FL  33620,  U.S.A. 
pwang@chumal.cas.usf.edu 


*U.S.  Army  Engineer 
Research  and  Development 
Center 

Coastal  and  Hydraulics 
Laboratory 

3909  Halls  Ferry  Road 
Vicksburg,  MS  39180,  U.S.A. 


ABSTRACT 


iittUUh. 


WANG,  P.  and  KRAUS,  N.C.,  2005.  Beach  profile  equilibrium  and  patterns  of  wave  decay  and  energy  dissipation 
across  the  surf  zone  elucidated  in  a  large-scale  laboratory  experiment.  Journal  of  Coastal  Research,  21(3),  522-534. 
West  Palm  Beach  (Florida),  ISSN  0749-0208. 

The  widely  accepted  assumption  that  the  equilibrium  beach  profile  in  the  surf  zone  corresponds  with  uniform  wave- 
energy  dissipation  per  unit  volume  is  directly  examined  in  six  cases  from  the  large-scale  SUPERTANK  laboratory 
experiment.  Under  irregular  waves,  the  pattern  of  wave-energy  dissipation  across  a  large  portion  of  the  surf  zone 
became  relatively  uniform  as  the  beach  profile  evolved  toward  equilibrium.  Rates  of  wave-energy  dissipation  across 
a  near-equilibrium  profile  calculated  from  wave  decay  in  the  surf  zone  support  the  prediction  derived  by  Dean  (1977). 
Substantially  different  equilibrium  beach-profile  shapes  and  wave-energy  dissipation  rates  and  patterns  were  gen¬ 
erated  for  regular  waves  as  compared  to  irregular  waves  of  similar  statistical  significant  wave  height  and  spectral 
peak  period.  Large  deviation  of  wave-energy  dissipation  from  the  equilibrium  rate  occurred  at  areas  on  the  beach 
profile  with  active  net  cross-shore  sediment  transport  and  substantial  sedimentation  and  erosion.  The  rate  of  wave- 
energy  dissipation  was  greater  at  the  main  breaker  line  and  in  the  swash  zone,  as  compared  to  middle  of  the  surf 
zone.  Based  on  analysis  of  the  SUPERTANK  data,  a  simple  equation  is  developed  for  predicting  the  height  of  irregular 
waves  in  the  surf  zone  on  an  equilibrium  profile.  The  decay  in  wave  height  is  proportional  to  the  water  depth  to  the 
one-half  power,  as  opposed  to  values  of  unity  or  greater  derived  previously  for  regular  waves. 
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INTRODUCTION 

The  concept  of  an  equilibrium  shape  of  the  beach  profile 
has  proved  fruitful  in  a  variety  of  applications  in  coastal  sci¬ 
ence  and  engineering,  such  as  beach  nourishment  design  and 
studies  concerned  with  morphological  evolution  in  the  near¬ 
shore.  Under  this  concept,  if  the  incident  waves  and  water 
level  remain  constant,  the  beach  profile  is  expected  to  evolve 
toward  a  stable  shape  for  which  the  net  cross-shore  sediment 
transport  rate  across  the  profile  approaches  zero.  Assump¬ 
tions  underlying  the  simplest  such  situation  involving  ana¬ 
lytical  solutions  are  alongshore  uniformity,  uniform  or  nearly 
uniform  sized  sediment,  absence  of  obstructions  such  as  hard 
bottom  on  the  beach  profile,  and  dominance  of  wave  action 
as  the  driving  process  shaping  the  profile.  These  assumptions 
can  be  relaxed  in  numerical  solutions,  such  as  that  of  Larson 
and  Kraus  (2000)  for  describing  profile  evolution  in  the  pres¬ 
ence  of  hard  bottom.  In  the  study  described  here,  we  inves¬ 
tigate  the  equilibrium  profile  shape  in  response  to  breaking 
waves,  i.e.,  the  profile  shape  in  the  surf  zone  on  sandy  beach- 
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es.  Inman  et  al.  (1993)  discussed  equilibrium  profile  shapes 
inside  and  outside  the  surf  zone,  with  the  offshore  bar  acting 
as  a  hinge  point. 

In  a  pioneering  study  of  the  beach  profile,  Bruun  (1954) 
noted  that  the  shape  of  the  nearshore  beach  profile  along  the 
coast  of  Denmark  followed  a  simple  power  law  of  depth  with 
distance  seaward  from  the  shoreline.  In  a  macro-scale  ap¬ 
proach  depending  on  general  characteristics  of  the  incident 
waves  and  representative  beach  material  grain  size,  Dean 
(1977)  derived  classes  of  beach  profile  shapes  that  depend  on 
choice  of  the  governing  mechanism  assumed  to  control  the 
transport  processes.  The  choice  that  best  corresponded  to 
field  data  for  Atlantic  Ocean  and  Gulf  of  Mexico  beaches 
yielded  a  shape  similar  to  that  found  by  Bruun  (1954)  in  a 
derivation  for  which  the  wave-energy-flux  dissipation  per 
unit  water  volume  is  uniform.  The  result  of  Dean  (1977)  has 
led  to  a  standard  equilibrium  beach  profile  shape  given  by: 

h  =  A/3  (1) 

where  h  is  the  still-water  depth,  x  is  the  horizontal  distance 
from  shoreline  ( h  =  0),  and  A  is  a  dimensional  scale  param¬ 
eter  having  units  of  m1/3  and  a  value  determined  mainly  by 
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sediment  grain  size  (Dean,  1977;  Moore,  1982;  Dean,  1991; 
Moutzouris,  1991). 

Many  studies  have  shown  that  Equation  (1)  or  a  similar 
form  represents  the  time-  and  locally-averaged  shape  of  the 
beach  profile  with  ample  sediment  supply  on  a  coast  exposed 
to  waves  and  uninfluenced  by  structures  and  hard  bottom 
(e.g.,  Moore,  1982;  Larson,  1991;  Work  and  Dean,  1991;  Pru- 
szak,  1993;  Larson  and  Kraus,  1994a;  Wang  and  Davis, 
1998,  1999;  Larson  et  al.,  1999;  Gonzalez  et  al.,  1999).  Lar¬ 
son  and  Kraus  (1989)  derived  an  analytical  form  for  the 
beach  profile  shape  at  equilibrium  under  the  realistic  wave 
decay  model  of  Dally  et  al.  (1985)  that  has  a  planar  shape 
at  the  foreshore  joining  to  the  form  of  Equation  (1)  further  in 
the  surf  zone.  Larson  et  al.  (1999)  derived  the  power-function 
beach  profile  (Eq.  1)  from  a  non-local  balance  of  time-aver¬ 
aged  offshore  and  onshore  sediment  transport.  These  and 
other  studies  give  confidence  that  Equation  (1)  both  describes 
nature  and  has  a  consistent  theoretical  basis.  However,  to  our 
knowledge,  no  equilibrium  mechanism  has  been  directly  test¬ 
ed  by  comparisons  to  laboratory  or  field  measurements  of 
wave  dissipation  and  grain  size.  Hereafter,  this  paper  consid¬ 
ers  wave-energy  dissipation  per  unit  volume  in  the  surf  zone 
and  its  relation  to  the  equilibrium  profile. 

In  the  Dean  (1977)  derivation  of  Equation  (1)  based  on  uni¬ 
form  dissipation  of  wave-energy  flux  per  unit  volume  across 
the  surf-zone  profile,  the  spilling  breaking  wave  assumption 
was  invoked,  meaning  the  height  H  of  broken  waves  in  the 
surf  zone  is  proportional  to  the  depth  h.  The  equilibrium 
wave-energy  dissipation  D,  across  the  surf  zone  is  then 
(Dean,  1977): 


Dt.  =  ^y2Pwg3/2As/2  (2) 

where  pw  is  the  density  of  water,  g  is  the  acceleration  of  grav¬ 
ity,  and  y  is  the  breaker  index  for  depth-limited  waves,  de¬ 
fined  as  wave  height  over  water  depth  ( H/h )  for  breaking  or 
broken  waves.  Kaminsky  and  Kraus  (1993)  found  that  the 
commonly  used  value  for  incipient  breaking  of  depth-limited 
individual  breaking  waves,  y  =  0.78,  is  supported  by  a  large 
laboratory  database  for  regular  (monochromatic)  waves.  They 
also  found  limits  for  y  in  the  more  than  400  individual  mea¬ 
surements  to  lie  within  0.60  and  1.59. 

For  a  train  of  waves  with  different  heights  (irregular 
waves),  at  a  given  location  in  the  surf  zone,  the  y  value  for 
significant  wave  height  Hs  is  also  found  to  be  close  to  0.78. 
However,  root-mean-square  wave  height  Hrms  is  appropriate 
for  calculation  of  mean  wave  energy  in  an  irregular  wave 
field.  At  a  given  location,  this  quantity  has  a  y  (rms)  value 
in  the  approximate  range  of  0.4  to  0.6  for  irregular  waves  in 
the  surf  zone  (e.g.,  Thornton  and  Guza,  1982,  1983;  Wang  et 
al.,  2002a,  2002b),  because  the  ratio  H/h  is  formed  by  the 
combination  of  breaking,  broken,  and  non-breaking  waves. 
Field  measurements  analyzed  by  Thornton  and  Guza  (1982, 
1983)  gave  a  value  of  yrms  in  the  range  of  0.42  to  0.44. 

Wave-energy-flux  dissipation  per  unit  water  volume  across 
shore  D(x)  can  be  calculated  from  measurements  of  wave 
height  and  water  depth  at  adjacent,  closely  spaced  wave 
gauges  as: 


h  ax  n 


d\  \pwgHrms^gh\ 


dx 


1  Pwg^  A(H2rmshV2) 
8  h,„„,  \x 


(3) 


where  E  is  the  wave  energy  per  unit  area,  Ce  =  \/gh  is  the 
wave  group  velocity  in  shallow  water,  and  hmid  is  the  water 
depth  at  the  mid-point  between  the  wave  gauges  used  here 
to  approximate  the  energy  dissipation  per  unit  volume. 

The  concept  of  beach-profile  equilibrium  and  the  corre¬ 
sponding  uniform  wave-energy  dissipation  underlie  a  number 
of  practical  cross-shore  sediment-transport  and  beach-profile 
evolution  models  (e.g.,  Kriebel  and  Dean,  1985;  Kriebel, 
1986;  Larson  and  Kraus,  1989;  Zheng  and  Dean,  1997).  In 
development  of  these  models,  it  is  assumed  that  wave-energy 
dissipation  rates  (Eq.  3)  will  deviate  substantially  from  the 
equilibrium  rate  (Eq.  2)  during  high-wave  energy  storms,  be¬ 
cause  the  equilibrium  beach  profile  and  associated  wave-en¬ 
ergy  dissipation  would  tend  to  represent  the  average  condi¬ 
tion  of  all  waves,  as  demonstrated  with  high-accuracy  field 
data  (Larson  and  Kraus,  1994a).  Therefore,  the  net  cross¬ 
shore  sediment  transport  rate  q  and  beach-profile  evolution 
depend  on  the  magnitude  of  deviation  from  an  equilibrium 
state  (e.g.,  Kriebel  and  Dean,  1985;  Kriebel,  1986;  Kraus 
and  Larson,  1988;  Larson  and  Kraus,  1989): 


q  =  K(D  -D,)»  (4) 

where  K  and  n  are  empirical  coefficients.  The  original  studies 
employing  this  form  took  n  =  1,  but  larger  values  such  as  n 
=  1.5  have  been  advocated  to  increase  the  time  rate  of  change 
at  longer  elapsed  time  (Zheng  and  Dean,  1997). 

Because  Equation  (2)  quantifies  the  equilibrium  rate  of  dis¬ 
sipation  per  unit  water  volume  in  both  theoretical  and  ap¬ 
plied  treatments  of  beach-profile  change,  it  is  of  central  in¬ 
terest  to  assess  its  validity  under  realistic  conditions.  The 
validation  should  answer  two  questions:  1)  Is  the  wave-en¬ 
ergy  dissipation  rate  uniform  across  the  equilibrium  surf- 
zone  profile  as  assumed,  and  2)  If  (1)  is  true,  how  well  does 
Equation  (2)  predict  the  magnitude  of  the  equilibrium  rate  of 
wave-energy  dissipation  per  unit  volume? 

The  SUPERTANK  large-scale  laboratory  data  collection 
project  (Kraus  et  al.,  1992;  Kraus  and  Smith,  1994;  Smith 
and  Kraus,  1995)  documented  beach-profile  response  to  ac¬ 
cretionary  waves,  erosional  waves,  and  different  boundary 
conditions,  together  with  high-density  measurements  of 
waves  and  currents  in  the  surf  zone.  The  SUPERTANK  data 
set  has  enabled  examination  of  a  wide  range  of  cross-shore 
processes  from,  for  example,  beach  profile  response  to  a  sea¬ 
wall  (McDougal  et  al.,  1996)  through  detailed  investigations 
of  net  cross-shore  sand  transport  (Larson  and  Kraus,  1994b) 
and  surf-zone  hydrodynamics  (Kriebel  and  Smith,  1994; 
Smith,  1994). 

For  the  present  study,  profile  survey  data  and  wave  mea¬ 
surements  across  the  surf  zone  made  at  SUPERTANK  are 
analyzed  to  test  the  assumptions  that  at  or  near  equilibrium 
of  the  beach  profile,  wave-energy-flux  dissipation  per  unit 
water  volume  is  uniform  across  the  surf  zone,  and  with  a 
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Table  1.  SUPERTANK  data  used  in  this  study. 


Run  ID 

Duration 

of  wave 

action 

(min) 

Hmo  before 
breaking 

(m) 

TP 

(s) 

7sw 

Comments 

ST-10-1-20 

20 

0.8 

3.0 

20 

initial  profile 

ST-10-1-270 

270 

0.8 

3.0 

20 

end  profile 

ST.30-1-9 

9 

0.4 

9.0 

20 

initial  profile 

ST.30-1-279 

279 

0.4 

9.0 

20 

end  profile 

ST-10-2-5 

5 

0.8 

4.5 

3.3 

initial  profile 

ST.10-2- 170 

170 

0.8 

4.5 

3.3 

end  profile 

ST.30-2-9 

9 

0.4 

8.0 

3.3 

initial  profile 

ST.30-2-209 

209 

0.4 

8.0 

3.3 

end  profile 

ST.GO-1-9 

9 

0.8 

3.0 

Mon 

initial  profile 

ST-G0-1-210 

210 

0.8 

3.0 

Mon 

end  profile 

ST-10-1-20 

20 

0.5 

8.0 

Mon 

initial  profile 

ST-10-1-559 

559 

0.5 

8.0 

Mon 

end  profile 

value  close  to  that  given  by  Equation  (2).  Six  cases  from  the 
SUPERTANK  dataset  were  selected  that  had  constant  waves 
run  for  at  least  3  hours.  The  cases  cover  considerably  differ¬ 
ent  wave  conditions.  The  objectives  of  this  study  are  to  ex¬ 
amine  (1)  patterns  of  surf-zone-profile  evolution  under  vari¬ 
ous  wave  conditions;  (2)  patterns  of  wave  decay  and  wave- 
energy  dissipation  associated  with  the  profile  change  and  in¬ 
cident  waves;  (3)  hypothesis  that  uniform  wave-energy-flux 
dissipation  per  unit  volume  occurs  over  an  equilibrium  pro¬ 
file;  and  (4)  the  accuracy  of  Equation  (2)  if  uniform  wave- 
energy  dissipation  is  achieved  on  an  near-equilibrium  beach 
profile.  In  addition,  a  new  equation  is  derived  that  describes 
decay  in  height  of  irregular  waves  on  an  equilibrium  profile. 

SUPERTANK  EXPERIMENTS 

SUPERTANK  was  designed  to  collect  data  without  scaling 
distortions  to  improve  understanding  and  predictive  capabil¬ 
ities  on  cross-shore  sediment  transport  processes  (Kraus  et 
al. ,  1992).  The  experiments  were  conducted  in  the  large  wave 
channel  at  the  0.  H.  Hinsdale  Wave  Research  Laboratory  at 
Oregon  State  University.  The  channel  is  104  m  long,  3.7  m 
wide,  and  4.6  m  deep,  into  which  a  76-m  long  sand  beach  was 
emplaced  for  the  SUPERTANK  experiments.  The  flap-type 
wave  generator  was  equipped  with  a  sensor  and  feedback 
mechanism  to  absorb  reflected  waves  at  the  peak  energy.  The 
beach  was  composed  of  approximately  600  m3  of  very-well 
sorted  fine  sand  with  a  median  diameter  of  0.22  mm  and  an 
average  fall  speed  of  3.3  cm/s.  Based  on  equations  developed 
by  Moore  (1982)  and  Dean  (1991),  the  A-value  in  Equation 
(1)  corresponding  to  this  grain  size  is  0.10  m1/3.  A  larger  A- 
value  of  0.17  m1/3  was  specified  in  the  initial  beach  construc¬ 
tion  to  create  a  steeper  profile  to  ensure  adequate  water 
depth  in  the  offshore  area  (Kraus  and  Smith,  1994).  In  design 
of  the  SUPERTANK  tests,  erosion  or  accretion  of  the  profile 
was  predicted  prior  to  running  the  waves  by  relations  given 
by  Kraus  et  al.  (1991)  to  arrange  instruments,  video  cameras, 
and  satisfy  requirements  for  participating  investigators. 

Six  cases  from  the  SUPERTANK  data  set  were  judged  suit¬ 
able  (with  sufficient  duration)  for  this  study  (Table  1).  In  the 
following  discussion  and  figures,  the  run  ID,  e.g.,  ST.10- 1- 
20,  indicates  the  following  information:  1)  the  first  five  char¬ 


acters,  ST.10,  denote  the  SUPERTANK  test  case,  (which  in¬ 
cludes  a  number  of  wave-run  segments  typically  with  10  to 
70  min  of  wave  action  each);  2)  the  next  two  characters,  -1, 
indicate  a  selected  block  in  the  ST.10  experiment,  which  is 
composed  of  several  bursts  or  time  intervals  with  the  same 
input  waves,  and  3)  the  characters  following  the  last  hyphen 
indicate  the  number  of  minutes  of  the  particular  wave  action 
(could  be  the  sum  of  several  wave-run  segments,  particularly 
for  the  end  profile).  Still-water  level  in  the  channel  was  con¬ 
stant  and  the  same  for  these  tests. 

Narrow-peaked  (swell  type)  irregular  waves  with  ysw  (spec¬ 
tral  width  parameter)  =  20  (Bouws  et  al.,  1985)  were  gen¬ 
erated  during  the  first  two  cases,  ST.10-1  and  ST.30-1. 
Steep  erosional  waves  with  a  significant  wave  height  Hs  of 
0.8  m  and  spectral  peak  wave  period  Tp  of  3  s,  were  generated 
in  ST.10-1.  The  initial  beach  profile  was  constructed  based 
on  Equation  (1)  with  an  A-value  of  0.17  m1'3.  The  end  profile 
resulted  from  a  total  of  270  min  of  wave  action.  Mild-steep- 
ness  accretionary  waves  with  a  significant  height  of  0.4  m 
and  peak  period  of  9  s  were  generated  in  ST.30-1.  The  initial 
profile  was  the  final  profile  of  a  previous  130-min  wave  run 
(Hs  =  0.5  m,  Tp  =  8  s,  ysw  =  3.3).  The  end  profile  of  ST.30- 
1  was  shaped  by  a  total  of  279  min  of  wave  action. 

Broad-peaked  (sea  type)  irregular  waves  with  ysw  =  3.3 
were  generated  during  the  cases,  ST.10-2  and  ST.30-2. 
Steep  erosional  waves  with  significant  height  of  0.8  m  and 
peak  period  of  4.5  s  were  generated  during  ST.10-2.  The  ini¬ 
tial  profile  was  the  final  profile  of  a  previous  280-min  narrow- 
peak  wave  run  ( Hh  =  0.8  m,  Tp  =  4.5  s,  y3W  =  20).  The  end 
profile  of  ST.10-2  was  shaped  by  a  total  of  170  min  of  wave 
action.  Gentle  accretive  waves,  with  significant  height  of  0.4 
m  and  peak  period  of  8  s  were  generated  in  ST.30-2.  The 
initial  profile  was  the  final  profile  of  a  previous  wave  run  of 
varying  wave  heights  and  periods.  The  end  profile  was 
shaped  by  a  total  of  209  min  of  wave  action. 

Regular  (monochromatic)  waves  were  generated  during  the 
last  two  cases  of  SUPERTANK,  ST.G0-1  and  ST.10-1.  Steep 
erosional  waves  with  height  of  0.8  m  and  period  of  3  s  were 
generated  in  ST.G0-1.  The  significant  wave  height  and  peak 
wave  period  are  similar  to  those  of  ST.10-1,  except  the  waves 
were  monochromatic.  Different  energy  levels  are  carried  by 
regular  and  irregular  waves  with  identical  significant  wave 
heights.  This  is  because  for  regular  waves,  significant  wave 
height  equals  root-mean-square  (rms)  wave  height,  which  is 
directly  proportional  to  wave  energy,  whereas  for  irregular 
waves,  significant  wave  height  equals  1.4  times  the  rms  wave 
height  for  Rayleigh-distributed  waves.  The  initial  profile  was 
the  final  profile  of  a  previous  wave  run  of  varying  wave 
heights  and  periods.  The  end  profile  of  ST.G0-1  resulted 
from  a  total  of  210  min  of  wave  action.  Mild-steepness  accre¬ 
tionary  waves  with  height  of  0.5  m  and  period  of  8  s  were 
generated  in  ST.10-1.  The  significant  wave  height  and  peak 
period  are  similar  to  those  of  ST.30.2,  except  that  the  waves 
were  monochromatic.  The  initial  profile  was  the  final  profile 
of  a  previous  wave  run  of  varying  wave  heights  and  periods. 
The  end  profile  of  ST.10-1  was  shaped  by  559  min  total  wave 
action. 

The  beach  profile  was  surveyed  across  the  middle  of  the 
3.66-m  wide  channel  at  the  end  of  each  wave-run  segment, 
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typically  10  to  70  min  long.  A  Geodimeter  T140  survey  sta¬ 
tion  was  positioned  outside  of  the  tank  to  track  by  an  infrared 
beam  the  position  of  the  prism  mounted  on  the  top  of  the  rod. 
Reproducibility  of  the  elevation  measurement  was  typically 
within  0.5  cm  except  at  the  shoreward  face  of  steep  break¬ 
point  bars,  where  the  variation  in  reproducibility  decreased 
to  a  nominal  2  cm  (Kraus  et  al. ,  1992). 

Waves  were  measured  at  26  cross-shore  locations  by  resis¬ 
tance  gauges,  spanning  from  the  swash  zone  to  approximate¬ 
ly  70  m  seaward  of  the  shoreline.  The  wave  gauges  were  sam¬ 
pled  at  16  Hz.  Both  spectral  analysis  and  time-series  zero- 
upcrossing  analysis  were  applied  to  calculate  wave  height 
and  period  (Kraus  and  Smith,  1994).  A  low-pass  filter  was 
applied  to  separate  incident-band  wave  motion  from  low-fre¬ 
quency  motions.  A  high-pass  time  series  was  thus  obtained 
by  subtracting  the  low-pass  time  series  from  the  preprocessed 
data  (Kraus  and  Smith,  1994,  pp. 62-66).  The  period  cutoff 
for  the  filter  was  set  to  twice  the  peak  period  of  the  incident 
waves.  This  data  filtering  has  considerable  influence  on  the 
calculation  of  wave  height  and  period  in  the  swash  zone, 
where  long-period  motion  was  significant  and  sometime  dom¬ 
inant.  The  wave  data  analyzed  in  this  paper  were  processed 
with  the  above  data  filtering  and,  therefore,  the  wave  height 
referenced  in  the  following  analysis  does  not  include  contri¬ 
butions  from  low-frequency  motions.  The  root-mean-square 
wave  height  and  mean  wave  period  calculated  from  the  zero- 
upcrossing  analysis  were  used  to  calculate  wave-energy  dis¬ 
sipation. 

The  rate  of  wave-energy-flux  dissipation  was  calculated 
from  the  measurements  at  two  adjacent  wave  gauges  based 
on  Equation  (3)  and  is  represented  (plotted  in  the  following 
figures)  at  the  mid-point  between  them.  Typically,  eight  wave 
gauges  were  located  in  the  swash  zone,  spaced  at  0.9-m  in¬ 
terval.  Two  wave  gauges  were  located  seaward  of  the  swash 
gauges,  spaced  at  1.8  m.  Sixteen  gauges  were  located  through 
and  beyond  the  surf  zone,  spaced  at  3.7-m  interval.  Given  the 
close  spacing  between  the  wave  gauges,  Equation  (3)  is  ex¬ 
pected  to  provide  reliable  estimates  of  wave-energy  dissipa¬ 
tion  rate  per  unit  water  volume  for  the  SUPERTANK  wave 
conditions,  except  possibly  in  areas  of  extreme  gradients  in 
wave  height,  e.g.,  in  the  vicinity  of  plunging  breakers. 

RESULTS 

In  the  six  cases  examined  and  nearly  all  the  SUPERTANK 
runs,  a  breakpoint  bar,  or  bar  in  the  area  of  incipient  depth- 
limited  wave  breaking,  developed  on  the  initial  monotonic 
power-function  beach  profile  as  described  by  Equation  (1). 
Little  change  occurred  in  the  offshore  region  outside  the  surf 
zone  in  the  cases  examined  here.  All  profiles  examined  in  this 
study  are  displayed  in  Figure  1.  The  surf-zone  portions  of  the 
profiles  are  not  directly  comparable  among  different  tests, 
particularly  in  the  vicinity  of  shoreline  where  the  beach 
might  have  been  reconfigured  for  different  objectives.  How¬ 
ever,  the  offshore  area  was  not  reconstructed.  The  following 
discussion  pertains  the  beach  profile  at  and  landward  of  the 
break-point  bar  (the  surf  zone). 


Figure  1.  All  beach  profiles  examined  in  this  study;  note  that  the  profile 
changed  little  seaward  of  the  bar. 


Beach-Profile  Evolution 

Test  ST.10-1  with  steep,  narrow-band  erosional  waves  was 
the  first  run  of  SUPERTANK.  The  initial  profile  was  the 
monotonic  power-function  profile  constructed  with  A  =  0.17 
m173.  Considerable  shoreline  recession  of  approximately  2.5  m 
and  development  of  a  30-cm  high  break-point  bar  occurred 
(Figure  2A).  The  erosion  in  the  vicinity  of  the  shoreline 
roughly  equaled  the  accumulation  in  the  bar  area,  while 
beach  elevation  remained  largely  unchanged  over  a  consid¬ 
erable  area  from  15  to  20  m  between  the  swash  zone  and  the 
main  breaker  line  over  the  bar  (Figure  2A).  This  morphologic 
change  indicates  that  the  sand  eroded  from  the  shoreline  was 
transported  offshore  across  an  area  in  the  mid-surf  zone  and 
deposited  near  the  main  breaker  line  to  form  the  break-point 
bar.  Moderate  erosion  in  the  trough  and  slight  accumulation 
directly  landward  of  it  also  occurred. 

The  initial  surf-zone  profile  represents  the  power-function 
profile  with  A  =  0.17  m1/3  well,  as  constructed.  The  construct¬ 
ed  profile  was  slightly  gentler  in  the  inner  surf  zone  and 
steeper  in  the  outer  surf  zone,  as  compared  to  A  =  0.17  m1/3 
profile.  After  270  min  of  wave  action,  the  end  surf-zone  pro¬ 
file  had  a  much  gentler  overall  slope  than  the  initial  slope 
because  of  erosion  in  the  vicinity  of  the  original  shoreline  and 
development  of  the  break-point  bar.  An  equilibrium  profile 
shape  with  A  =  0.17  m1/3  over-predicts  the  water  depth  across 
most  of  the  surf  zone,  except  at  the  bar  trough.  The  value  A 
=  0.10  m1/3,  which  corresponds  to  the  SUPERTANK  grain 
size  of  0.22  mm  according  to  Moore  (1982)  and  Dean  (1991), 
predicts  the  profile  shape  reasonably  well  across  most  of  the 
surf  zone,  except  at  the  trough,  where  water  depth  is  under¬ 
predicted  considerably. 

For  the  remainder  of  the  three  irregular  wave  cases  (Figure 
2B  and  Figure  3),  the  magnitude  of  the  beach-profile  change 
was  much  smaller  as  compared  to  the  initial  changes  (Figure 
2A),  especially  during  the  accretionary  wave  cases  with  lower 
waves  (Figure  2B  and  Figure  3B).  The  bar  became  slightly 
wider  and  flatter,  as  compared  to  the  initial  profile  after  the 
accretionary  wave  actions.  Considerable  berm  erosion  above 
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Figure  2.  Initial  and  end  beach  profiles  for  ST_10-1  and  ST-30-1,  and 
the  power-function  profiles. 


mean  water  level  occurred  during  the  steep,  broadband  (sea 
type)  erosional  waves  (Figure  3A).  The  eroded  sand  was  de¬ 
posited  in  a  relatively  thin  layer  across  the  mid-surf  zone 
extending  to  the  bottom  of  the  trough.  The  average  elevation 
from  8  to  25m  across  shore  gained  2.4  cm  during  the  170-min 
wave  action.  The  elevation  of  the  bar  remained  nearly  stable 
at  the  beginning  and  end  of  the  wave  action,  but  moved  ap¬ 
proximately  0.5  m  offshore. 

Similar  to  ST.10-1,  the  surf-zone  profile  can  be  represent¬ 
ed  more  accurately  with  A  =  0.10  m1/3  (except  at  the  trough), 
whereas  A  =  0.17  m1/3  over-predicts  water  depth  across  most 
of  the  surf  zone  except  at  the  trough,  especially  for  the  ac¬ 
cretive  wave  cases  (Figures  2B  and  3B).  In  other  words,  the 
initially  constructed  steep  profile  was  out  of  equilibrium  and 
evolved  toward  a  gentler  profile,  as  predicted  by  the  DEAN 
(1977)  profile  model.  The  A  of  0.10  m173  gives  an  equilibrium 
profile  shape  matching  the  end  profile  well. 

Considerable  change  between  the  initial  and  end  profiles 
occurred  during  the  regular  wave  cases  (Figure  4A  and  B). 
During  the  erosional  case  with  higher  waves  (Figure  4A),  the 
bar  moved  onshore  and  became  considerably  higher  after  210 
min  of  wave  action.  Moderate  erosion  near  the  shoreline  also 
occurred.  No  profile  change  occurred  seaward  of  the  bar  out¬ 
side  the  surf  zone,  whereas  considerable  changes  occurred 
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Figure  3.  Initial  and  end  beach  profiles  for  ST_10-2  and  ST.30-2,  and 
the  power-function  profiles. 


throughout  the  surf  zone.  The  surf-zone  profile  was  not  as 
smooth,  with  several  secondary  bars  (or  large  ripples)  appar¬ 
ent,  as  compared  to  the  profile  during  the  irregular  cases. 
The  power-function  model  (Eq.  1)  with  A  =  0.17  m1/3  over- 
predicts  the  water  depth  across  most  of  the  surf  zone,  where¬ 
as  the  lower  A  =  0.10  m173  under-predicts  the  water  depth  at 
most  locations.  The  overall  surf-zone  profile  is  substantially 
different  from  that  under  irregular  waves,  although  the  reg¬ 
ular  wave  case  applied  similar  statistical  significant  wave 
height  and  peak  wave  period  run  in  the  irregular  case. 

For  the  accretionary  regular  waves  with  smaller  height 
(Figure  4B),  substantial  profile  changes  occurred  after  559 
min  of  wave  action.  The  bar  moved  onshore  substantially 
while  preserving  height.  Slight  advance  of  the  shoreline  oc¬ 
curred.  The  inner  surf  zone  landward  of  the  break-point  bar 
gained  a  substantial  amount  of  sand.  Large  bedforms,  or  sec¬ 
ondary  bars,  developed  on  the  seaward  slope  of  the  bar.  The 
equilibrium  profile  shape  with  A  =  0.10  m173  slightly  over- 
predicts  most  of  the  profile  except  at  the  trough,  whereas  A 
=  0.17  m173  substantially  over-predicts  the  water  depth  across 
nearly  the  entire  profile. 

The  overall  beach  profile  shape  under  regular  waves  was 
different  as  compared  to  the  profile  shape  corresponding  to 
the  irregular  waves  of  similar  statistical  significant  wave 
height  and  peak  period.  The  erosional  regular  wave  produced 
a  steeper  profile,  whereas  the  accretionary  regular  wave  re¬ 
sulted  in  a  gentler  profile.  Also,  secondary  bars  or  large  rip¬ 
ples  developed  in  the  mid-surf  zone  under  regular  waves,  in 
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Figure  4.  Initial  and  end  beach  profiles  for  ST.GO-l  and  STJO-1,  and 
the  power-function  profiles. 


contrast  to  the  smooth  profile  under  irregular  waves.  The 
beach  profile  in  the  surf  zone  remained  reasonably  similar 
under  irregular  waves  with  different  spectral  characteristics, 
as  well  as  statistical  wave  height  and  wave  period.  Constancy 
of  profile  shape  supports  present  understanding  that  the  A- 
value  defining  the  equilibrium  beach-profile  shape  is  not  sig¬ 
nificantly  influenced  by  wave  height,  period,  and  spectral 
characteristics,  leaving  sediment  grain  size  as  the  dominant 
factor  controlling  the  beach-profile  shape,  as  found  in  nu¬ 
merous  previous  studies  using  field  data  (e.g.,  Dean,  1977; 
Moore,  1982;  Moutzouris,  1991). 

Patterns  of  Wave  Decay  and  Breaker  Index  in  the  Surf 
Zone 

In  this  study,  the  main  breaker  line  for  both  regular  and 
irregular  wave  cases  was  defined  to  be  at  the  location  land¬ 
ward  of  which  a  significantly  steeper  rate  of  wave-height  de¬ 
cay  was  measured.  This  criterion  was  based  on  comprehen¬ 
sion  that  substantial  wave-energy  loss  or  wave-height  de¬ 
crease  should  follow  major  wave  breaking.  The  rate  of  wave 
decay  at  the  breaker  line  may  also  provide  information  on 
the  intensity  of  wave  breaking.  For  example,  WANG  et  al. 
(2002a)  measured  a  much  greater  rate  of  wave-height  de¬ 
crease  following  plunging  breaking  than  that  following  spill¬ 
ing  breaking  in  a  mid-scale  laboratory  facility. 

Different  patterns  of  wave  decay  were  measured  in  ST.10- 
1  at  the  beginning  and  the  end  of  the  wave  run  (Figure  5A). 
At  the  breaker-point  bar,  a  steeper  decay  in  wave  height  oc¬ 


Figure  5.  Cross-shore  distributions  of  wave-height  decay  and  breaker 
index  over  the  initial  and  end  beach  profile  for  ST_10-1  and  ST_30-1. 


curred  at  the  end  of  the  wave  run  than  at  the  beginning.  This 
evolution  in  pattern  of  wave  decay  follows  the  evolution  of 
the  breakpoint  bar  and  change  in  its  seaward  slope.  Smith 
and  Kraus  (1991)  showed  that  waves  breaking  on  bars  have 
different  characteristics  than  those  breaking  on  uniformly 
sloping  beaches.  In  the  vicinity  of  the  shoreline,  a  gentler 
decay  was  measured  at  the  end,  probably  influenced  by  the 
concave  beach-profile  shape  there.  Wave  conditions  seaward 
of  the  surf  zone  remained  nearly  identical  during  the  270  min 
wave  action. 

The  breaker  index  (BI)  yrms  =  HrmJh  for  irregular  waves  in 
the  surf  zone  exhibited  an  increasing  trend  in  value  with  dis¬ 
tance  from  the  breakpoint  landward,  from  0.37  at  the  breaker 
line  to  1.05  in  the  swash  zone  at  the  beginning  of  the  wave 
run  (Figure  5A).  At  the  end,  a  greater  breaker-index  value  of 
0.60  was  measured  at  the  breaker  line,  followed  by  a  decrease 
landward  to  0.45.  Thereafter,  a  steady  increasing  trend  to 
approximately  0.88  in  the  swash  zone  was  measured,  mean¬ 
ing  that  virtually  all  waves  were  breaking  or  had  broken  in 
that  swash  zone. 

The  patterns  of  wave  decay  at  the  beginning  and  end  of  the 
wave  action  for  the  other  three  irregular  wave  cases  were 
nearly  identical  (Figure  5B;  Figures  6A  and  B),  especially  for 
the  two  accretionary  cases  with  relatively  small  wave  heights 
(Figures  5B  and  6B).  Similar  wave  patterns  are  consistent 
with  the  stable  beach  profiles.  The  slight  broadening  and  flat- 
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Figure  6.  Cross-shore  distributions  of  wave-height  decay  and  breaker 
index  over  the  initial  and  end  beach  profile  for  ST.10-2  and  ST.30-2. 


tening  of  the  bar  did  not  result  in  any  substantial  change  in 
wave  pattern  over  the  bar.  Breaking  of  the  infrequent  larger 
waves  probably  caused  the  slight  flattening  over  the  bar, 
which  did  not  induce  significant  change  in  the  overall  wave 
height  and  wave-height  decay  over  the  bar. 

The  breaker  index  yrma  increased  from  0.20  over  the  sea¬ 
ward  slope  of  the  bar  to  0.75  in  the  swash  zone  (Figures  5B 
and  6B).  The  decrease  of  Hrma  and  Haig  over  the  bar  was  not 
as  steep  as  that  in  the  higher  wave  cases  (Figures  5A  and 
6A),  indicating  that  wave  breaking  over  the  bar  was  not  as 
significant.  Steep  wave  decay,  indicating  major  wave  break¬ 
ing,  occurred  close  to  the  shoreline. 

For  the  erosional  wave  case  ST.10-2,  the  main  breaker  line 
as  indicated  by  a  steep  decay  of  wave  height  was  located  over 
the  bar  crest  (Figure  6A).  At  the  beginning  of  the  wave  run, 
the  wave  decay  in  the  vicinity  of  the  shoreline  was  slightly 
gentler  than  that  at  the  end  of  the  run.  The  slight  seaward 
movement  of  the  bar  during  the  170  min  wave  action  did  not 
result  in  any  distinguishable  changes  in  the  patterns  of  wave 
decay  (Figure  6A). 

The  breaker  index  yrma  exhibited  similar  cross-shore  distri¬ 
bution  pattern  as  that  of  ST.10-1,  with  a  peak  at  the  bar, 
followed  by  a  moderate  decrease,  and  then  a  landward  in¬ 
creasing  trend  (Figure  6A).  Over  the  bar,  the  breaker  index 


Figure  7.  Cross-shore  distribution  of  -yrmI. 


equaled  roughly  0.65,  followed  by  a  decrease  to  0.44  over  the 
trough,  and  then  increased  to  1.10  in  the  swash  zone. 

Overall,  the  breaker  index  yrma  showed  a  landward  increas¬ 
ing  trend  for  all  the  four  irregular  wave  cases.  An  increased 
yrma  value  over  the  breaker-point  bar  was  produced  by  the 
combination  of  an  increased  wave  height  due  to  shoaling  and 
the  shallow  water  or  the  bar.  Figure  7  compiles  all  yrma  values 
for  the  four  cases.  The  horizontal  axis  represents  the  dimen¬ 
sionless  distance  to  the  still-water  shoreline,  with  x  being  the 
cross-shore  distance  to  shoreline,  and  xb  is  the  width  of  the 
surf  zone.  This  landward  increasing  trend  of  the  breaker  in¬ 
dex  across  the  surf  zone  is  described  reasonably  well  with  a 
power  function  as: 

/xV0-36  (xYV3 

7ml8  =  0.416- 

\xb)  \xbJ 

(I  \~2/9 

k)  (5) 

where  the  term  with  the  power  -0.36  represents  the  least- 
square  power-function  fit  of  the  measured  yrma  with  the  cor¬ 
relation  coefficient  R2  equals  0.68.  The  right  side  of  Equation 
(5)  was  obtained  by  substitution  with  Equation  (1),  and  the 
representative  value  of  yrmaJl  fits  well  in  the  range  of  0.42- 
0.44  reported  by  Thornton  and  Guza  (1982,  1983)  for  field 
measurements. 

Combining  Equations  (5)  and  (1),  a  simple  equation  for  pre¬ 
dicting  the  decay  of  Hrma  across  the  surf  zone  of  an  equilib¬ 
rium  beach  profile  is  obtained  as 

Hrma  =  y  rms_bAx1,3x%3  =  yrma_bx%3\/Ah 

=  yrms.b^Kh  (6) 

where  yrmaJ,  is  the  value  of  yrma  at  the  point  of  incipient  depth- 
limited  wave  breaking.  The  function  at  the  second  and  third 
lines  of  the  equation  is  obtained  by  substituting  the  power- 
function  beach  profile  (Eq.  1)  into  the  first  line.  Care  should 
be  exercised  in  determining  yrmaJ>  over  a  barred  profile  and 
might  best  be  defined  past  the  plunge  area  of  incipient  break¬ 
ing  waves.  For  comparison,  THORNTON  and  GUZA  (1983) 
obtained  the  result  Hrma  ~  h9'10  for  irregular  wave  breaking 
on  a  plane  sloping  bottom  as  h  — »  0  in  the  surf  zone. 


Journal  of  Coastal  Research,  Vol.  21,  No.  3,  2005 


Beach  Profile  Equilibrium  and  Wave-Energy  Dissipation 


529 


Figure  8.  Measured  and  predicted  (Eq.  6)  wave  heights  across  the  surf  zone,  irregular  wave  cases. 


Measured  Hrms  values  were  compared  with  the  predicted 
ones  from  Equation  (6).  The  surf-zone  width  xb  was  measured 
from  the  still-water  shoreline  to  the  breaker  line  where  rel¬ 
atively  steep  wave-energy  dissipation  occurred.  The  breaking 
depth  hb  is  determined  based  on  Equation  (1)  from  the  mea¬ 
sured  xb.  The  depth  hb  can  also  be  measured  from  the  beach 
profile.  However,  measured  hb  is  not  directly  used  here  be¬ 
cause  it  is  sensitive  to  the  relief  of  the  bar.  The  best-fit  yrmsJb 
value  of  0.42  (Eq.  5)  was  used  in  the  calculation.  More  data 
are  necessary  to  improve  the  reliability  and  accuracy  of  a 
general  ■yrm3_b. 

The  Hrms  was  predicted  reasonably  well  for  the  two  ero- 
sional  cases  with  higher  waves  (Figure  8:  two  left  panels). 
The  overall  trend  of  wave  decay  was  predicted  well  although 
the  wave  height  was  under-predicted  at  most  locations  except 
in  the  swash  zone.  A  10%  increase  in  yrmsJ(  value  would  im¬ 
prove  the  match  with  the  measured  data.  The  wave  height 
over  the  breaker-point  bar  was  substantially  under-predicted, 
and  a  discrepancy  is  expected  because  Equation  (6)  pertains 
to  a  monotonic  beach  profile. 

For  one  of  the  two  accretionary  cases,  ST.30-1  with  nar¬ 
row-peak  wave  and  long  period  of  9  s,  wave-height  increases 
were  measured  at  two  locations,  a  short  distance  landward 
of  the  bar  and  directly  seaward  of  the  swash  zone.  These  are 
likely  due  to  wave  reformation  after  breaking  on  the  bar  and 
shoaling  before  further  breaking  in  the  swash  zone.  Although 


the  complex  variations  of  wave  height  were  not  predicted  by 
Equation  (6),  the  general  trend  was  reproduced  reasonably 
well  (Figure  8:  upper  right  panel).  Equation  (6)  is  developed 
based  on  the  analysis  of  the  monotonic  equilibrium  beach  pro¬ 
file  as  given  by  Equation  (1)  and  cannot  reproduce  compli¬ 
cated  changes  in  wave  height  over  a  multiple  bar  profile. 

Wave-height  decay  was  reproduced  well  for  the  other  ran¬ 
dom  wave  accretionary  case,  ST.30-2,  except  for  over-predic¬ 
tion  at  the  barred  breaker-line  (Figure  8:  bottom  right  panel). 
Wave  decay  at  the  breaker  line  for  the  accretionary  wave  case 
with  lower  wave  height  was  much  gentler  as  compared  to  the 
erosional  cases.  The  minor  peak  consisting  of  a  reduced  rate 
of  wave  decay  followed  by  an  increased  rate  in  the  mid-surf 
zone  was  not  predicted. 

The  patterns  of  wave  decay  differed  at  the  beginning  and 
the  end  of  the  erosional  regular  wave  case,  ST.G0-1  (Figure 
9A).  At  the  beginning,  a  steady  decrease  in  wave  height  with 
approach  to  shore  was  measured,  whereas  much  steeper 
wave  decay  was  measured  over  the  bar  after  the  210  min 
wave  action.  The  large  wave-height  decay  from  approximate¬ 
ly  1.1  m  to  0.4  m  over  the  bar  apparently  resulted  from  the 
uniform  breaking  of  regular  waves  over  a  narrow  zone,  in 
contrast  to  wider  breaker  zones  of  irregular  waves.  The  wave 
height  remained  almost  constant  across  most  of  the  surf  zone, 
followed  by  a  steep  decrease  at  the  shoreline.  Some  wave- 
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Figure  9.  Cross-shore  distributions  of  wave-height  decay  and  breaker 
index  over  the  initial  and  end  beach  profile  for  ST.G0-1  and  STJ0-1.  The 
straight  line  represents  =  0.62  (solid  line)  and  y  =  0.78  (dashed  line). 


height  variations  were  measured  seaward  of  the  surf  zone, 
probably  resulting  from  formation  of  standing  waves. 

At  the  beginning  of  the  wave  run,  the  breaker  index  re¬ 
mained  relatively  uniform  across  shore,  typically  ranging 
from  0.53  to  0.87.  At  the  end  of  the  wave  action,  the  value  of 
the  breaker  index  was  nearly  1.20  at  the  breaker  line,  fol¬ 
lowed  by  a  sharp  drop  to  0.24  over  the  trough  and  then  a 
steady  increase  to  approximately  1.04  in  the  swash  zone.  The 
standard  7  =  0.78  found  in  regular  wave  studies  (Kaminsky 
and  Kraus,  1993)  yielded  over-predictions  across  most  of  the 
surf  zone,  except  at  the  breaker  line  and  in  the  swash  zone 
(Figure  9A).  Time-series  wave  analysis  yielded  an  average 
HrmJHsig  ratio  of  0.80  in  the  surf  zone  for  both  of  the  two 
regular  wave  cases.  This  value  is  greater  than  the  0.71  ob¬ 
tained  from  the  Rayleigh-distribution  assumption.  Outside 
the  surf  zone,  the  HrmJHsig  ratio  averaged  0.96,  which  is  close 
to  unity  as  expected  for  regular  waves.  A  sharp  decrease  in 
the  HrmJHsig  ratio  from  nearly  unity  to  approximately  0.8  oc¬ 
curred  in  the  vicinity  of  the  breaker  line,  or  the  narrow  zone 
of  initial  wave  breaking.  The  HrmJHsig  ratio  of  0.8  would  yield 
a  yrms  value  of  0.62  given  the  standard  7Big  =  0.78.  The  value 
7„  =  0.62,  also  plotted  in  Figure  9A,  yielded  an  improved 
match  with  the  measurements,  with  an  average  yrma  value  of 
0.66. 

Substantially  different  wave-decay  patterns,  as  compared 
to  the  irregular  wave  cases,  were  also  observed  in  the  accre¬ 


tionary  regular  wave  case,  ST.I0-1.  The  main  breaker  line 
as  indicated  by  the  steep  wave-height  decay  moved  onshore, 
consistent  with  the  substantial  onshore  movement  of  the  bar 
(Figure  9B).  Somewhat  regulated  variations  of  the  wave 
height,  perhaps  due  to  formation  of  partial  standing  waves, 
occurred  both  inside  and  outside  the  surf  zone.  The  breaker 
index  varied  from  0.36  to  over  1.41  across  the  surf  zone.  The 
landward  increasing  trend  was  not  as  steady  as  in  the  other 
cases.  A  breaker  index  greater  than  1.0  was  measured  over 
the  bar  crest,  i.e.,  at  the  main  breaker  line,  both  at  the  be¬ 
ginning  and  at  the  end  of  the  wave  run.  The  standard  7  = 
0.78  matched  the  measured  values  (with  an  average  of  0.74) 
better  than  the  7rma  =  0.62.  The  reason  for  this  is  not  clear. 
Substantially  different  patterns  and  trends  of  wave-height 
decay  were  produced  during  regular-wave  experiments,  as 
compared  to  the  results  with  irregular  waves  of  similar  sta¬ 
tistical  offshore  significant  wave  height  and  peak  period.  The 
landward  increasing  trend  of  yrms  during  the  regular  wave 
cases  was  not  as  clear  as  that  during  the  irregular  wave  cas¬ 
es.  Equation  (5)  developed  for  irregular  waves,  does  not  de¬ 
scribe  the  regular  wave  situation.  A  much  greater  value  of 
yrms,  exceeding  1.0,  was  measured  at  the  breaker  line  over 
the  bar  for  the  regular  wave  cases.  As  found  from  several 
previous  regular  wave  studies  (Dally  et  al.,  1985;  Smith  and 
Kraus,  1987;  Smith  and  Kraus,  1991),  the  wave-height-decay 
curve  tends  to  be  concave  (curved  or  rounded  downward). 
This  is  in  contrast  to  the  present  irregular  wave-decay  curve, 
which  has  the  sense  of  convex  upward  (somewhat  rounded 
upward).  The  difference  indicates  that  results  obtained  from 
simpler  (in  terms  of  wave  generation)  regular  wave  studies 
might  not  be  directly  applicable  to  irregular  waves,  similar 
to  the  situation  of  beach-profile  change  as  discussed  above. 
Wave  irregularity  should  be  considered  in  the  modeling  of 
wave-height  decay  and  beach-profile  evolution  in  the  surf 
zone. 

Patterns  of  Wave-Energy  Dissipation  Per  Unit  Volume 

The  key  assumption  of  the  Dean  (1977)  equilibrium  theory 
is  that  wave-energy  dissipation  per  unit  volume  is  uniform  if 
the  beach  profile  is  in  equilibrium.  In  the  following,  the  as¬ 
sumption  is  examined  with  the  SUPERTANK  data. 

Significantly  different  patterns  of  wave-energy  dissipation 
per  unit  volume  were  measured  at  the  beginning  and  end  of 
ST.10-1  (Figure  10A).  At  the  beginning  of  the  run,  a  sub¬ 
stantial  landward  increasing  trend  was  measured  over  the 
initial  monotonic  power-function  profile.  The  wave-energy 
dissipation  rate  per  unit  volume  increased  from  nearly  zero 
just  seaward  of  the  surf  zone  to  over  400  Nm/m3/s  in  the 
swash  zone.  After  270  min  of  wave  action,  the  rate  of  wave- 
energy  dissipation  per  unit  volume  became  much  more  uni¬ 
form,  ranging  from  54  to  155  Nm/m3/s  with  an  average  of  107 
Nm/m3/s  and  a  standard  deviation  of  43  Nm/m3/s.  The  equi¬ 
librium  energy  dissipation  rate  calculated  from  Equation  (2) 
( A  =  0.1  m1/3;  7rms  =  0.56)  was  65  Nm/m3/s,  approximately 
39%  lower  than  the  average  value  of  the  measured  dissipa¬ 
tion  rate.  The  value  0.56  for  yrms  was  obtained  by  dividing 
the  standard  7  value  of  0.78  (for  significant  wave  height)  by 
V2,  on  the  assumption  of  a  Rayleigh  distribution  in  wave 
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Figure  10.  Wave-energy  dissipation  rate  per  unit  volume  over  the  initial 
and  end  beach  profiles,  and  the  computed  equilibrium  dissipation  rate  for 
ST.10-1  and  ST.30-1. 


height.  At  the  end  of  the  wave  action,  a  high  rate  of  wave- 
energy  dissipation  was  measured  over  the  bar  crest  at  the 
main  breaker  line,  followed  by  a  sharp  drop  over  the  trough 
and  a  gradual  increase  toward  the  swash  zone. 

Although  the  measured  yrmB  values  were  not  uniform  across 
the  surf  zone  (as  discussed  in  the  previous  section)  and  on 
average  were  greater  than  0.56,  they  were  not  used  in  the 
Equation  (2)  calculation,  because  the  equation  was  developed 
under  the  assumption  of  a  constant  y  under  spilling  breakers. 
Predictions  from  Equation  (2)  based  on  the  simple  equilibri¬ 
um  profile  theory  matched  the  measured  values  within  50% 
for  this  case. 

Based  on  the  overall  beach-profile  shape  and  the  relatively 
minor  profile  changes  in  a  series  of  following  runs,  it  is  rea¬ 
sonable  to  believe  that  the  end  profile  after  the  270  min  of 
wave  action  should  be  close  to  the  equilibrium  shape.  Large 
deviations  from  both  measured  and  predicted  equilibrium  en¬ 
ergy  dissipation  occurred  near  the  shoreline  and  at  the  bar 
crest  at  the  beginning  of  the  ST.10-1  wave  run.  These  devi¬ 
ations  corresponded  to  substantial  beach-profile  change  and, 
therefore,  net  cross-shore  sediment  transport.  The  results 
qualitatively  support  the  transport  equations  (e.g.,  Eq.  4)  in 
the  profile-change  models  of  Kriebel  and  Dean  (1985)  and 
Larson  and  Kraus  (1989). 

For  the  two  irregular  accretionary  wave  cases,  the  observed 
wave-energy  dissipation  was  nearly  the  same  at  the  begin¬ 


ning  and  end  of  wave  action  (Figures  10B  and  11B).  This 
result  implies  nearly  identical  wave-decay  pattern  and  beach- 
profile  shape,  indicating  that  the  beach  was  already  in  near¬ 
equilibrium  state  at  the  beginning  of  the  run.  Due  to  the  rel¬ 
atively  small  wave  height,  wave  breaking  over  the  bar  crest 
may  not  be  as  significant  as  the  high-wave  cases,  as  also  in¬ 
dicated  by  the  small  wave-energy  dissipation  rate  at  the  bar 
crest.  Major  wave  breaking  occurred  in  a  narrow  zone  near 
the  shoreline,  as  indicated  by  increased  dissipation  of  wave 
energy  there.  Overall,  measured  rates  of  wave-energy  dissi¬ 
pation  were  moderately  greater  than  the  predictions  from 
Equation  (2)  near  the  shoreline.  In  contrast,  across  the  rest 
of  the  surf  zone,  the  predicted  rate  was  greater  than  the  mea¬ 
sured  rates,  which  were  only  slightly  above  zero.  Some  wave 
regeneration  was  evident,  particularly  for  ST.30-1  (the  nar¬ 
row-peaked  swell  type  accretive  wave  case),  as  indicated  by 
energy  production  (negative  value  of  wave-energy  dissipa¬ 
tion).  The  3.7-m  spacing  may  also  limit  accuracy  of  calcula¬ 
tion  of  energy  dissipation  where  gradients  in  wave  height  are 
complicated  or  where  waves  are  generated  at  the  plunge 
point.  On  average,  the  wave-energy  dissipation  rate  for  the 
two  accretionary  cases  was  46  Nm/m3/s  for  the  narrow-peak 
swell-type  waves  and  44  Nm/m3/s  for  the  broad-peak  sea-type 
waves.  The  predicted  value  of  65  Nm/m3/s  (Eq.  2)  is  44% 
greater  than  the  average  measured  value,  although  the  mea¬ 
sured  values  showed  considerable  variation  across  shore  with 
a  slight  increasing  trend  toward  the  shoreline. 

For  the  broadband  erosional  wave  case,  the  patterns  of 
wave-energy  dissipation  per  unit  volume  were  similar  at  the 
beginning  and  end  of  the  wave  action,  ranging  from  53  to  173 
Nm/m3/s  (Figure  11A).  Little  change  in  the  beach  profile  oc¬ 
curred,  indicating  that  it  was  probably  near  equilibrium  at 
the  beginning  of  the  run.  A  great  rate  of  wave-energy  dissi¬ 
pation  per  unit  volume  was  measured  over  the  bar  crest  at 
the  main  breaker  line,  followed  by  a  sharp  decrease  and  a 
slight  increase  toward  the  swash  zone.  A  sharp  increase  of 
the  dissipation  rate  occurred  in  the  swash  zone.  The  wave- 
energy  dissipation  rate  per  unit  volume  across  most  of  the 
mid-surf  zone  was  relatively  uniform  and  matched  well  with 
the  value  obtained  from  Equation  (2).  The  two  peaks  in  dis¬ 
sipation  at  the  bar  crest  and  in  the  swash  zone  were  more 
than  twice  the  predicted  value  from  Equation  (2). 

Different  patterns  of  wave-energy  dissipation  rate  per  unit 
volume  were  measured  during  the  two  regular  wave  cases, 
as  compared  to  the  irregular  wave  cases.  Substantially  dif¬ 
ferent  patterns  were  also  measured  at  the  beginning  and  the 
end  of  the  wave  run,  consistent  with  the  different  beach-pro¬ 
file  shapes  and  wave-height  decay.  At  the  beginning  of  the 
regular  erosive  wave  case,  the  wave-energy  dissipation  rate 
per  unit  volume  ranged  from  —43  to  402  Nm/m3/s  with  no 
apparent  trend  of  changing.  A  negative  wave-energy  dissi¬ 
pation  rate  indicates  an  increase  in  wave  height.  This  is  ap¬ 
parent  for  the  accretionary  regular  wave  case.  At  the  end  of 
the  210  min  wave  action,  a  high  dissipation  rate  of  nearly 
624  Nm/m3/s  was  obtained  over  the  bar  crest  at  the  main 
breaker  line,  followed  by  a  sharp  decrease  to  24  Nm/m3/s  just 
landward  of  the  trough  (Figure  12A).  The  wave-energy  dis¬ 
sipation  rate  remained  close  to  zero  across  most  of  the  mid¬ 
surf  zone,  and  then  increased  to  nearly  292  Nm/m3/s  in  the 
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Figure  11.  Wave-energy  dissipation  rate  per  unit  volume  over  the  initial 
and  end  beach  profiles,  and  the  computed  equilibrium  dissipation  rate  for 
ST.10-2  and  ST_30-2. 
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Figure  12.  Wave-energy  dissipation  rate  per  unit  volume  over  the  initial 
and  end  beach  profiles,  and  the  computed  equilibrium  dissipation  rate  for 
ST.G0-1  and  ST J0-1. 


swash  zone.  The  pattern  of  wave-energy  dissipation  rate  per 
unit  volume  deviated  greatly  from  uniformity  and  differed 
significantly  from  the  predicted  value  from  Equation  (2).  De¬ 
viation  from  uniformity  of  wave-energy  dissipation  rate  was 
much  greater  at  the  end  of  the  wave  run  than  at  the  begin¬ 
ning,  influenced  by  the  intense  wave  breaking  over  a  narrow 
zone  over  the  bar. 

For  the  regular  wave  accretionary  case,  the  energy  dissi¬ 
pation  rate  per  unit  volume  ranged  from  —204  to  637  Nm/ 
m3/s  at  the  beginning  of  the  wave  action,  which  is  substan¬ 
tially  different  from  the  —296  to  751  Nm/m3/s  at  the  end  (Fig¬ 
ure  12B).  The  dissipation  patterns  were  far  from  being  uni¬ 
form  and  also  differ  substantially  from  the  predicted  value  of 
65  Nm/m3/s. 

The  modulation  of  standing  waves  might  have  significant 
influence  on  the  pattern  of  wave-energy  dissipation  for  the 
regular  wave  cases,  especially  for  the  long-period  accretion¬ 
ary  wave  case,  as  indicated  by  the  large  negative  dissipation 
rate  in  the  surf  zone.  The  extremely  high  peak  of  wave-en¬ 
ergy  dissipation  at  the  main  breaker  line  in  both  of  the  reg¬ 
ular  wave  cases  is  apparently  related  to  the  concentrated 
wave  breaking  over  the  bar. 

The  more  gradual  breaking  of  irregular  as  compared  to  reg¬ 
ular  waves  tended  to  form  a  smoother  beach  profile  and 
wave-decay  pattern,  whereas  the  localized  breaking  of  regu¬ 


lar  waves  produced  more  irregularity  in  the  beach  profile  and 
shape  of  wave  decay.  Given  the  steadiness  of  the  monochro¬ 
matic  waves,  concentrated  wave  breaking  and  energy  dissi¬ 
pation  should  be  expected,  especially  for  plunging  type  break¬ 
ing.  The  deviation  of  HrmJHaig  from  unity  in  the  surf  zone  for 
regular  waves  indicates  a  certain  amount  of  wave  irregularity 
is  generated  in  the  breaking  process.  No  substantial  differ¬ 
ences  were  found  between  the  narrow-banded  swell-type 
waves  and  broad-banded  sea-type  waves  in  the  data  exam¬ 
ined  in  this  study. 

CONCLUSIONS 

In  the  six  SUPERTANK  cases  analyzed,  a  nearshore  bar 
developed  and  was  maintained  over  an  initially  monotonic 
power-function  profile.  Except  in  the  direct  vicinity  of  the 
break-point  bar,  the  power-function  model  of  Dean  (1977)  de¬ 
scribed  the  equilibrium  beach  profile  reasonably  well.  For  ir¬ 
regular  waves,  the  assumption  that  an  equilibrium  profile 
shape  develops  under  uniform  wave-energy-flux  dissipation 
is  supported  by  the  SUPERTANK  data  across  a  large  portion 
of  the  surf  zone.  Relatively  greater  rates  of  wave-energy  dis¬ 
sipation  were  measured  at  the  main  breaker  line  and  in  the 
swash  zone.  The  equilibrium  rate  of  wave-energy  dissipation 
was  predicted  mostly  within  50%  by  the  simple  model  derived 
by  Dean  (1977)  with  a  spilling  breaking  assumption.  Signif- 


Journal  of  Coastal  Research,  Vol.  21,  No.  3,  2005 


Beach  Profile  Equilibrium  and  Wave-Energy  Dissipation 


533 


icant  net  cross-shore  sediment  transport  and  beach-profile 
change  occurred  at  locations  where  the  rate  of  wave-energy 
dissipation  deviated  substantially  from  the  equilibrium  rate. 

For  irregular  waves,  in  the  surf  zone  the  root-mean-square 
wave  height  on  an  equilibrium  beach  profile  was  found  to  be 
proportional  to  the  square  root  of  local  water  depth,  except 
over  the  break-point  bar.  The  more  gradual  wave-height  de¬ 
cay  contrasts  to  the  common  spilling  wave  breaker  assump¬ 
tion  for  which  the  wave  height  is  proportional  to  local  water 
depth  in  the  surf  zone.  The  simple  model,  requiring  the  input 
of  the  equilibrium  A  value  and  the  bounds  of  the  surf  zone, 
reproduced  the  measured  wave  decay  reasonably  well  under 
different  large-scale  laboratory  wave  and  beach  conditions. 

Both  beach-profile  evolution  and  wave-energy  dissipation 
pattern  differ  substantially  under  regular  waves  as  compared 
to  those  under  irregular  waves  of  comparable  statistical  sig¬ 
nificant  wave  height  and  peak  period.  Different  beach  profile 
shapes  were  developed  under  regular  waves,  particularly  in 
the  vicinity  of  the  break-point  bar.  For  the  regular  wave  cas¬ 
es,  uniformity  of  wave-energy  dissipation  did  not  improve  be¬ 
tween  the  start  and  end  of  wave  actions.  These  results  indi¬ 
cate  studies  with  regular  waves  may  not  represent  all  sedi¬ 
ment-transport  processes  as  associated  with  irregular  waves. 
A  key  difference  is  that  wave-energy  dissipation  for  regular 
waves  tends  to  be  concentrated  in  a  narrow  zone  of  breaking. 
Knowledge  of  cross-shore  sediment  transport  gained  from 
regular-wave  experiments  will  carry  uncertainties  if  applied 
to  irregular  waves  without  appropriate  interpretation  or 
treatment  of  the  regular  waves,  such  as  through  a  Monte 
Carlo  approach. 

For  the  irregular  wave  cases  analyzed,  there  were  relative¬ 
ly  small  differences  between  breaker  index  values,  wave  de¬ 
cay,  and  energy-dissipation  rate  per  unit  volume  for  the 
broad  and  narrow-banded  waves.  This  result,  based  on  lim¬ 
ited  laboratory  data  but  at  large  scale,  indicates  that  the  in¬ 
fluences  of  spectral  characteristics  of  irregular  waves  on  the 
beach-profile  evolution  are  not  significant  for  typical  waves 
with  single-peaked  spectra. 
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